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Abstract 

We show that the existence of exceptional polynomials leads to the presence of non-trivial super- 
symmetry. The existence of these polynomials reveals several distinct isospectral potentials for the 
Schrodinger equation. All Schrodinger equations having Laguerre and Jacobi polynomials as their solu- 
tions, have non-trivial supersymmetric partners with corresponding exceptional polynomials as solutions. 

1 Introduction 

The celebrated theorem of Bochner [1] states that if an infinite sequence of polynomials Pn (x) — y„ (n = 
0, 1, 2 • • • oo) satisfies a second order eigenvalue equation of the form, 

P{x)yn + (l{x)y'n + 'r{x)yn = XnVn, (1) 

then q{x) and r{x) are polynomials of order 2, 1 and respectively. This theorem has been extended 
to (/-difference equations by Askey [5]. The polynomials y„ form a complete set with respect to a positive 
measure. Till recently it was thought that only the classical orthogonal polynomial systems (OPS) like, the 
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Hermite, Laguerre and the Jacobi satisfy eq.(IT]). Note that yo is a constant in all these systems. It was shown 
[3J S] that it is possible to construct an OPS, which starts with (n = 1, 2 • • •) by a suitable modification 
of the weight function which forms a complete set. The newly constructed OPS was shown to satisfy a 
Sturm-Liouville equation of the form given in eq.(IT]), where r{x) is not a constant. Only the Laguerre and 
the Jacobi polynomials were shown to allow such extensions, which are known as exceptional polynomials. 
One can construct the exceptional Xi-Laguerre Polynomials C'^{x), k > 0, using the Gram-Schmidt proce- 
dure from the sequence O |4] , 

vi=x + k + l; Vi = {x + k)\ i > 2, (2) 



using the weight function 



X e 



defined in tlie interval x G (0, oo) and the scalar product 

poo 

{f,g)k^ / dxWk{x)f{x)g{x). (4) 







The weight function for the normal Laguerre polynomial Wk{x) — x'^'e^^, is multiplied by suitable factors 
such that one obtains a new Wk (x) such that one can construct the new OPS which does not have the zero 
degree polynomial. The exceptional Xi-Laguerre differential equation is 

Tk{y) = Ay, (5) 

where A = n — 1 and 

Tk{y) = -xy" + (^1^) [(fc + X -I- 1)?/ - y\. (6) 

For more details we refer the reader to [31 H] and the references therein. 

Similarly the exceptional Xi- Jacobi polynomials V^'^\x), for a ^ P and real, are obtained from the 
sequence 

ui = X — c, Ui = {x — by, * > 2, (7) 

where 

a = — a); 6 = ^-i-^, and c — b+—. (8) 
The scalar product is defined in the range [—1, 1] with the weight function 

^".^(^) • (9) 

As in the case of exceptional Laguerre case, the weight function of this new OPS is a rational extension of 
the classical Jacobi weight function, Wa^pix) = (1 — x)"{l + x)'^ . 
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These exceptional Xi-Jacobi polynomials obey the eigenvalue equation 

{x^ - l)y" + 2a (^^) [{x - c)y' -y] = Ay, (10) 

where A = (n — l)(a + /3 + n). 

It is well known that the solutions of the Schrodinger eigenvalue problems, with exactly solvable (ES) 
potentials involve classical OPS. The form of the complete solution is ^/W{x)yn, where W{x) is the weight 
function corresponding to the OPS, y„(n = 0, 1, 2 . . .), such that W {x)yn{x)ym{x)dx = Smn- The weight 
function is determined from the boundary condition of the given potential. By adding suitable terms to 
the potential, we can modify the weight function, so that we obtain exceptional polynomials as solutions to 
the Schrodinger equation (SE). Wc make a crucial observation that a quantum mechanical problem which 
admits classical Laguerre/ Jacobi polynomials as a solutions of the SE, after a modification of the potential 
will admit exceptional Laguerre/ Jacobi polynomials as solutions having same eigenvalues, with the ground 
state missing in one of the Hamiltonians. 



2 Exceptional Polynomials in Quantum Mechanics 

Theorem : By adding an extra term Ve{x) to the Laguerre/ Jacobi differential equation and demanding the 
solutions to be g{x) = and g{x) = (^3^ for the Laguerre and Jacobi respectively, where f{x) satisfies 

Xi- exceptional differential equation for the Laguerre and Jacobi respectively, Ve{x,m) can be determined 
uniquely. 

Proof : Let g{x) = L^{x) satisfy the Laguerre differential equation 

cP d 
x-^g{x) + (m + 1 - x)—g{x) + (A - m)g{x) = 0, (11) 

where A is an integer. By adding an extra term Ve{x,m) to the Laguerre differential equation: 

x-^9{x) + (m + 1 - x)—g{x) + (A + Ve{x,m) - m)g{x) = 0, (12) 

and setting g{x) = j^^^, where f{x) satisfies the Xi exceptional Laguerre differential equation 

- xf'ix) + (l^) [{m + x + l)f'{x) - fix)] = {n- l)f{x), (13) 
determines Ve{x,m) to be 

1 2m 

Ve{x,m) = -—---——^. (14) 

(x + m) [x + m)^ 



If g{x) = (^J^^y with f{x) satisfying the X — j exceptional differential equation 



= {n-j)f{x), (15) 
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one obtains Ve{x, m) to be 

Ve (x, m) = ■— — - - — — — — (16) 

[x + m) [x + my 

for the general case. Similarly for the Jacobi polynomials, g{z) = Pn ^ 1/2. s+a 1/2)^^^^ satisfying the 
Jacobi differential equation 

(1 - z^)g"{z) + + + 2)z].g'(z) + n(n + a + ^ + l)g{z) = 0, (17) 

adding an extra potential 14 (z) to it gives 

(1 - z^)g"{z) + [P-a-{a + P + 2)z]g'(z) + n{n + Ve{z) + a + /3 + 1)5(2) = 0. (18) 

Setting g{z) — (^^jy and demanding that f{z) satisfy the Xi exceptional Jacobi differential equation 



(z2 _ l)/"(^) + 2a (^) [(^ - c)nz) - f{z)] = A/(z), (19) 



where A = (n — l)(a + /3 + n). This determines Ve{z) to be 



2 26 



This ends the proof of our theorem. 

Example(l): For 3D oscillator case in natural units the radial equation is solved by 



i?(0=e^e-ii^(0 (21) 
where K{(^) satisfies Laguerre differential Equation 

^^^(^) + + ^ - ^)^^(^) + (A - 3 - 2l)K{i) = 0, (22) 

with f ~ . where A = 27t, + 3, n = l + 2n' and n' ~ 0, 1, 2.... and A = 2E, which gives E = n+ 1. Identifying 
with the eq. PT|) one has m ^ I + ^ and X = n' + I + ^. Then one gets 



1 2(g^ 



= TT-^ITK - TT^WKT- (23) 



Thus one obtains the modified 3D oscillator potential, studied in |^, [2, by adding V°'''^{£,,1) to the 3D 
oscillator potential which was. A quantum Hamilton- Jacobi analysis of the modified oscillator potential has 
been done in [llj . Example(2) As is well known, the radial equation for the Coulomb potential in natural 
units is 



(f , , 2d , , 



A 1 + 



R{r) = 0. (24) 



The solution to R{r) is of the form 

R{r) ^r^e-iK{r), (25) 
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where K{r) satisfies Laguerre differential equation 

r^K{r) + {2{l + 1) - r)^K{r) + (X - I - l)K{r) = 0. (26) 
dr-^ ar 

Thus, K{r) = L^i^iir) comparing this with eq. ([TT|) one gets X = n' + I, m ~ 21 + 1 and A = ^=55- Now, 
by replacing A by X + Ve{r,l) then the extra term to eq. and eg. ([25]) takes the form 

1 2(2/ + 1 

(r + 27TT) " (r + 2/ + 1) 

Example(3) For the Morse potential [T(T, the radial equation is solved by 



R{y) = y^--e--^yK{y), (28) 
where K{y) satisfies Laguerre differential equation 



y—^iy) + (2(5 - n) + 1 - y)—Kiy) + (A - (s - n))if (j/) = 0. (29) 



d 

' dy^ \»/ ' \ \ / ■ dy 
Here X(?/) = L^'*"^"(j/) with y = ^e"^, and s = A/a and the eigenvalues, X^ - {A - na)^. Then the 
exceptional part is given by 

s) = , . - (30) 

(y + s — n) (y + s — n)^ 

The solution of the Schrodinger equation with extra potential Ve{x) added for bound state problems form a 



complete set. Note that solutions are of the form ^yW{x)C^ this gives a proof for the completeness of the 
exceptional polynomials. 

It should be noted that all these potentials admitting the Laguerre differential equation (fTTI) can be 
brought to the same form by a point canonical transformation. They only differ in the constant values of A 
and ni. For example, in natural units = — E, by making a change of variable y = x^, one gets the 
Coulomb potential V^C'°"' — ^L— = _^ _ without altering the centrifugal term. Similarly, starting from 
the Morse potential, V ^ A^ + exp(— 2Q;a;) — 2B{A + a/2) exp(— ax), by change of variable r = exp(aa;) 
goes over to the Coulomb potential. 

A similar analysis goes through for the Jacobi polynomials. Consider the Scarf potential [S1[B]: 

V{x) = -A^ + {A^ +B'^ - Aa) sec^{ax) - B{2A - a) tan(ax), (31) 

whose solution is 

i;iz) = (1 - zY'-^^/\l + z)(^+^)/2p^^-^-i/2.-+A-i/2)(^)^ (32) 

where z — sin(ax), s = A/a and A = B/a. Here, pO^^^~'^/^^^+^ ^^^\z) = y are the classical Jacobi 
polynomials which satisfy the Jacobi differential equation 

(1 - z^)y" + [l3~a~{a + l3 + 2)z]y' + n{n + a + (3 + l)y = 0. (33) 
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Adding an extra potential Ve{z) and demanding the solution to be of the form, 



with p^^~-^~^^^''^~^^~^/'^'> (^z) satisfying the Xi-Jacobi differential equation eq.((TU]), determines Ve{z) to be 



> ~ 2A-1-2BZ (2A - 1 - 2Bz)2 ' ^"^^^ 
The Scarf potential in eq. pip has the Jacobi polynomials as solutions and all the potentials having Jacobi 
polynomials as solutions can be mapped to the Scarf potential by a suitable change of variable. 



3 Supersymmetric Quantum Mechanics 

The discovery of exceptional polynomials led Quesne to guess the superpotential W{x), demanding isospec- 
trality [B]. Supersymmetric quantum mechanics assumes isospectral Hamiltonians. In the previous section, 
we have constructed the two sets of isospectral Hamiltonians, of which, one Hamiltonian denoted by iJ^ 
with potential has the Laguerre/ Jacobi solution and another Hamiltonian denoted by with potential 
has exceptional Laguerre/ Jacobi solution. Therefore, it is natural ask, "can the two potentials thus 
constructed be partner potentials? Can one construct a superpotential, mahcalW{x), from which these two 
can be constructed?" Here V~ cannot be obtained from through the shape invariant arguments. Below, 
we answer these questions in the affirmative. 

First, we briefly review conventional SUSY and refer the reader to P| and [TU' for more details. In 
supersymmetry, the superpotential W{x) is deflned in terms of the intertwining operators A and as 

A=^+W{x), = -^+Wix). (36) 

ax ax 

This allows one to define a pair of factorized Hamiltonians as 

H+ = A^A = --^ + V+ix)-E, (37) 

H- = AA^ = —^ + V-ix)~E, (38) 

where E is the factorization energy. 

The partner potentials (x) are related to W{x) by 

V^{x) ^W^{x)Ty^'ix) + E, (39) 

where prime denotes differentiation with respect to x. The eas. (j37p and p8L imply 

H+A^ = A^H-, Ah+ = H-A. (40) 
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From the above one can see that the operators A and A'' act as a intertwining operator. These operators 
ahows one to go from wave function ItJj^) to the wave function \ip~) and vice versa. Our aim is to construct 
the operator A for the isospectral Hamihonians given in the previous section. We would also like to add that 
the conventional methods for obtaining the superpotential from the ground state eigenfunction [71 [12] does 
not work here. In our case, the wave functions for both the Hamiltonians are known and one can construct 
the operator A from the following relation 

using the operator O [3J [31 [^ , which connects the ordinary Laguerre polynomials to the exceptional Laguerre 
polynomials 

dLt-'{^)^Ct+^{x), (42) 

where (!? = (§+ ||)(y§ ~ oo) — oo. The superpotential, W{x), can be obtained by replacing ^ in A in terms 
of O. The superpotential yV(x) determined to be 

4 -j^' (^^) 

We obtain the intertwining operators defined in eq. (j36p . which take one from IV'jj^) to |V'i7) smd vice- versa, 
which involve classical and exceptional Laguerre polynomials. Thus, the form of these intertwining operators 
is universal. Hence, we prove the existence of exceptional polynomials leads to the presence of non-trivial 
supersymmetry. 

It should be noted that to find the isospectral potentials, it is not necessary that the Hamiltonians be 
of the form A^A and AA^ . It turns out that in the case of 3D oscillator, since the Hamiltonian is positive 
definite, the supersymmetry machinery goes through and the partner Hamiltonians are of the form A^ A 
and AA^ and we recover the results of Quesne [6 , for appropriate identification of the variables. By taking 
Vi+{x) = ia;2 + _ E and V-{x) ^ V,Ii(x) + Ve{x) then one obtains 

2w'{r) = vr' = v+ -V- = + (44) 



X 



{x^ + fc)2 (.t2 + k) ' 



From this one gets the superpotential W'{x) for the 3D oscillator. However, in the Coulomb problem one 
has negative eigenvalues. This makes the Hamiltonian negative definite. This can be overcome by taking 
—A^A to be B^B, which is positive semi-definite. Or one can always add a zero point energy as done in 
conventional SUSY, i.e Hq + e = e + E' , where E' is the conventional energy of the bound state Hydrogen 
atom which can be written as A = ^ , that makes Ho + e positive semi-dcfinitc without changing the wave 
function. The superpotential of the Coulomb problem is obtained by mapping it to the oscillator problem, 

with the following change of variable r — x^. One gets 2W'{x) to be 

1/°^^ /II 91c 

2W'{x) = ^ = _ ' (45) 
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Similar method will work for ES potentials which have classical Jacobi polynomials as solutions. Here 
the operator, which allows us to go from classical to exceptional Xi-Jacobi polynomials is given as 

0,P:-'^^+\x) = 2{fi -a){P + n)V:-f,, (46) 

where Oj — [a + /3 — (/? — q;)^] ((1 + + /? + l) + (/3 — Q!)(1 + a;). Taking the wave functions given in 
Eqs. ([5^ and (|34p as il>^{x) and V' (2^) respectively, we can construct W{x), which allows us to construct 
the intertwining operators whose form is universal. We recover Quesne's results for the Scarf potential 
constructed in ref [B]. 

4 Conclusion 

We have proved a theorem that whenever the Schrodinger equation is associated with the Laguerre/ Jacobi 
polynomials as solutions, one can also construct the exceptional Laguerre/ Jacobi polynomials as solutions 
for the Schrodinger equation by adding an extra term Ve{x) to the original potential. The form of the 
potential is universal and it depends on the extra terms, l/{x -\- k) for the Laguerre and — 6) for the 
Jacobi polynomials. We have also proved this theorem for the general exceptional Laguerre/ Jacobi 
polynomials. We have shown that these two set of Hamiltonians satisfy new kind of supersymmetry called 
nontrivial supersymmetry. 
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